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fT^ ' Abstract. We further discuss the Unipolar Inductor Model (UIM) coupled to GW 



o 



C/3 



emission (Dall'Osso et al. 2005) and compare it to observed properties of the two 



O ■ candidate ultrashort period binaries RX J0806+15 and RX J1914+24 . 

' We consider the measured orbital periods, period derivatives and inferred X-ray 

^ ' luminosities of these two sources and find constraints on system parameters in 

5— ( ' order for the model to account for them. We find that these properties point to 



the two sources being in different regimes of the UIM, with the requirement of 
low magnetic moment primaries (~ 10^" G cm"^) for both. 

Given this weak magnetization, RX J0806+15 has a sufficiently low luminosity 



M ' that it can be interpreted as having a primary spin almost synchronous to and 

just slightly slower than the orbital motion. Its measured orbital spin-up is only 
slightly affected by spin-orbit coupling and is mostly due to GW emission. 
RX J1914+24, on the other hand, is too bright in X-rays and has too slow an 
orbital spin-up for the same regime to apply. We suggest that this binary system 
may be emitting GWs at a significantly higher rate than implied by its measured 
tOo — 6 X 10~^^ rad s~^. The latter is explained, in this framework, by the primary 
spin being slightly faster than the orbital motion (q < 1.1). In this case, the 
associated spin-orbit coupling transfers to the orbit a significant amount of angular 
momentum, thus partially balancing that lost to GW emission. 
All expectations can be tested in the near future to confirm the viability of the 
model. 
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1. Introduction 

The nature of the two X-ray sources RX J1914+24 and RX J0806-I-15 has been long de- 
bated since their discoveries and the origin of their emission is still quite unclear. Several 
models have been proposed so far, none of which seems able to explain all of their ob- 
served properties without facing major difficulties (see Cropper et al. 2003 for a review 
on models and observations). 

The long-term timing analyses of Strohmayer (2004,2005), Hakala et al. (2003, 2004) 
and Israel et al. (2004), the most recent ones based on ^ 10 yrs of monitoring, have 
unambiguously shown that the single periods found in the power-spectra of both objects 
(respectively ~ 9.5 and 5.3 min) are steadily decreasing over time. The very nature of 
such periodicities is debated, evidence being presently in favour of its interpretation as 
the orbital period of a binary system. Indeed, no other significant periods are found in 
power spectra which is quite unexpected if, for example, the observed periods reflected 
the spins of accreting compact objects. 

If the observed periods are true orbital modulations, then it is extremely difficult to acco- 
modate Roche-lobe filling non-degenerate secondary stars orbiting a compact companion 
in such close binary systems, because of the extremely small masses implied. For this 
reason, double-degenerate binary systems represent the most likely option. 
As orbital periods are observed to decrease, accretion models in a double degenerate sys- 
tem are apparently ruled out because (conservative) accretion in these cases should lead 
to a period increase. Therefore, although with the above caveats, the double-degenerate 
scenario with no mass exchange seems at least more suited than others to explain the 
timing properties of these sources. 

In this context the Unipolar Inductor Model (UIM), originally proposed for the Jupiter-Io 
system (Goldreich & Lynden-Bell 1969), has been suggested as an alternative explana- 
tion to the observational properties of the above mentioned sources (Wu et al. 2002). In 
this version of the model two white dwarfs form a close binary system with a moderately 
magnetized primary and a secondary that, for practical purposes, is non-magnetic. If the 
primary spin (wi) is not perfectly synchronous with the orbital motion Wo (the secondary 
is tidally locked to the orbit) an e.m.f. is induced across the secondary star, idealized as 
a perfect conductor crossing the primary's magnetic field lines. In the presence of free 
charges, the e.m.f. induces a current flow between the two stars, along the sides of the 
flux tube connecting them. Currents in the atmosphere of the primary cross field lines 
at a depth where the electrical conductivity {a) becomes isotropical; at this point, they 
close the circuit returning back to the secondary (Goldreich and Lynden-Bell 1969, Li, 
Ferrario and Wickramashige 1998, Wu et al. 2002). 

Resistive dissipation of currents and the associated heating, that occur essentially in the 
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primary atmosphere, cause the observed soft X-ray emission. The source of the emission 
is uhimately represented by the mechanical energy of the relative motion between the 
primary spin and the orbit. Therefore, dissipation consumes this energy and brings the 
primary and the orbit into synchronous rotation. 

As shown by Dall'Osso et al. (2005, paper I from here on), perfect synchronization is 
never reached in the UIM applied to ultrashort period DDBs. The emission of GW con- 
tinuously injects energy in the electric circuit when lui < uJq and can thus sustain a 
permanent, slight asynchronism of the primary spin. We called this the system's steady- 
state asynchronism: binaries in this regime, i.e. with sufficiently short orbital periods, 
will have a primary spin just slightly slower than the orbital motion, the small asynchro- 
nism causing a continuous flow of electric current which is ultimately powered by GWs. 
In this work, we first address and clarify some general issues discussed in paper I. Then, 
we apply our model to the two candidate ultrashort period binaries RX J0806+15 and RX 
J1914-I-24. We obtain simple relations between a source luminosity, its measured rate of 
orbital spin-up, component masses and primary magnetic moment, that give constraints 
on system parameters for the model to be applicable to both sources. 
Detailed and systematic calculations and evolutionary implications are beyond the scope 
of this work and will be addressed elsewhere. 



2. The Unipolar Inductor Model: summary 

We summarize here briefly the features of the UIM, referring to Wu et al. (2002) and 
paper I for a more detailed description. The key ingredient of the model is that the 
primary spin (wi) differs from the orbital motion (wq); ot = uji/uIq is the asynchronism 
parameter. 

In an asynchronous system with orbital separation a the secondary star moves across 
the primary magnetic field lines with a relative velocity v = a{u)o — ui). The electric 
field induced through the secondary is. E = ^ ^ The associated e.m.f. is $ = 2R2E, 
R2 being the secondary's radius. The induced flow of currents is resistively dissipated 
essentially in the primary atmosphere, causing significant local heating and thus power- 
ing the soft X-ray emission. The Lorentz torque on cross-fleld currents in the primary 
atmosphere and in the secondary causes spin-orbit coupling, redistributing angular mo- 
mentum between the primary spin and the orbital motion. However, given the significant 
GW-emission expected from two white dwarfs with ultrashort orbital period, orbital an- 
gular momentum is continuously lost from the system, which in turn requires spin-orbit 
coupling to continue: this keeps the electric circuit active at all times. 
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Concerning the specific assumptions of this work, we consider binary systems hosting 
two degenerate white dwarfs, with the following mass-radius relation (Nauenberg 1972): 

M ( M 



= 0.0112 



(1) 



,1.4337 V 1-433, 
where M is expressed in solar masses. 

In paper I we have shown that tidal synchronization of the primary component is not 
expected to be efficient on the orbital evolutionary timescale, while it should well be effi- 
cient for a lower mass companion, in most cases of interest. The system orbital evolution 
is thus given by (eq. E4 in Appendix E of Wu et al. 2002): 

where Egr = -(32/5)G/c^[g/(l + g)]^Mfa^u;^ is the energy loss rate through GW emis- 
sion (Landau & Lifshitz 1951), the second term in parentheses represents the contribution 
of spin-coupling and W is the electric current dissipation rate (source luminosity). The 
function g{ijJo) < (see appendix A for its definition) represents two thirds of the orbital 
(kinetic plus gravitational) energy of the binary system, times a coefficient just slightly 
smaller than unity. The latter accounts for the additional tiny amount of orbital angular 
momentum that is continuously lost to keep the secondary into synchronous rotation. 
GWs give by definition a positive contribution to lOo while the contribution of electric 
coupling (W > by definition) depends on the sign of (1 — q) . Given the measured orbital 
spin-up of both sources under study, it must be assumed for them that, if a > 1, then 
\Egr\ > \W/{l-a)\. 

The quantity W can be expressed as (paper I): 



W 



3? ~ V c i [GMi(l-Fg)]iV6 {H/Ad)j{e) 



= k iJo' (1 - af (3) 

where 5R is the system effective resistance (see Wu et al. 2002), jii the primary's magnetic 
moment, a the height-averaged WD atmospheric conductivity, H the atmospheric depth 
at which currents cross magnetic field lines and return back to the secondary and Ad is 
the thickness of the arc-like cross section of the ciirrent layer at the primary atmosphere. 
The last equality defines the system's constant k. 
Finally, we recall the evolution equation for a (see paper I) : 
a W w„ 



a a{l — a)Iiu}g cUo 
where = a/a. 



(4) 



3. UIM coupled to GW emission: energy budget and the general solution 

In paper I we discussed the conditions that cause the current flow and, thus, the soft 
X-ray emission, never to be quenched by synchronization of the primary star. In the time- 
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independent, linear approximation, the evolution equation for a is solved by a simple 
exponential function, whose asymptotic value ctoo < 1 confirms this expectation. 
Numerical integration of the full evolution equation confirms that, as the orbital period 
changes over time, a does indeed evolve asymptotically towards unity, never reaching it 
exactly. We showed that, in the general solution, the value < 1 is attained, which 
conserves the total mechanical energy resevoir (Euim) of the circuit 

EuiM ^{l/2)Ii{ijl-cjl). (5) 

Therefore, by requiring that the time derivative of Euim is zero, we obtain (paper I) 

«oo(l - «oo) = 11/3 - (6) 

This expression for is implicit, as it depends on the source's actual value of 0;^, a 
function of a itself^. Despite this, the above formula is of great practical use as it relates 
the model parameter a™ to measured quantities (cjo, Uo). We will return on this in §131 
Given the importance of the steady-state solution, it is necessary to specify here its 
physical meaning and definition, starting from a discussion of the overall energy budget 
of the system. 

3.1. Asynchronism parameter and efficiency of spin-orbit coupling 

As a preliminary step, we summarize the relevant properties of spin-orbit coupling in the 
UIM. First recall the expression of the Lorentz torque Nl — W/[ujo{^ — «)] (paper I). 
This same torque Nl acts - with opposite sign - on both the primary star and the orbit. 
Therefore, angular momentum is conserved by spin-orbit coupling: all that is subtracted 
from one component is transferred to the other one. 
The rate of work done by Nl on the orbit is 
W 

while that done on the primary spin is: 



E^'^ = -NluJo = (7) 



E,p^n ^ Nl^, ^ = -aEi^'^'K (8) 

The sum Espin + -E^"''''' = ~W, which means that not all of the energy extracted from 
one component is transferred to the other one. The energy lost to ohmic dissipation is 
the energetic cost of spin-orbit coupling. 
The following picture thus results from the above formulae: 

— When a > 1 the circuit is powered at the expenses of the primary's spin energy. Only 
a fraction of this energy is transferred to the orbit, the rest being lost to ohmic 
dissipation. 



^ The problem liEis a non-linear nature. 
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— When a < 1 the circuit is powered at the expenses of the orbital energy. Now a 
measures directly the fraction of energy that is transferred to the primary spin, while 
(1 — a) is that lost to ohmic dissipation^. 

Therefore, the parameter a can be seen as a measure of the efficiency with which spin- 
orbit coupling transfers energy between the orbit and the primary spin: the more asyn- 
chronous a system is, the less efficiently energy is transferred, most of it being dissipated 
as heat. 



3.2. Stationary state and the general solution 

We can now re-discuss and clarify the energetic argument of paper I, based on which 
we obtained the expression of . The steady-state regime is characterized by having a 
constant electric energy reservoir in the circuit, for a given lUo- The degree of asynchronism 
is thus set by the condition that the system uses up just the power fed to the circuit by 
GWs and spin-orbit coupling. It is important to note, in this context, that both ohmic 
dissipation and the energy required to change Wo and/or loi must be taken into account. 
Therefore, the condition Euim = of paper I would not be wholly consistent, because 
only LUi and lUq appear in Euim and no account is taken of the luminosity W. We then 
re-state the problem of steady-state in the following terms: since the circuit dissipates 
energy at a rate W, it must recharge at the same rate as the binary evolves, if its energy 
reservoir is to be conserved. The proper steady-state condition is thus 

Euim = ±W , (9) 

where the sign on the right-hand side depends on the asynchronism parameter. The plus 
sign holds when Ejjim > and, thus, a > 1, while the minus sign is relevant to the case 
with a < 1 {Euim < 0). 

In paper I (Appendix B) we have shown that the case with aoo > 1 corresponds to a 
condition very hardly realized in WD+WD systems, because it can hold only for very 
close systems. So close, indeed, that the lighter component would most likely fill its Roche- 
lobe. Even if it did not, however, the condition also requires lUo < 0, which means the 
system must widen; the two components will rapidly be too far apart for the condition to 
hold. Therefore, proper steady-state exists in WD+WD systems only with Euim ~ —W. 
Inserting the expressions for uji and W in condition Q we obtain the expression for the 
steady-state asynchronsim. We call it simply aoo, to distinguish it from the (slightly 
different) expression of a™ of paper I. 

Our former assumption Euim — amounts ultimately to neglecting W compared to 
the energy terms contained in Euim ■ Note that the luminosity W is produced by ohmic 
^ Note that, for a = 0, Espin = but Espin > 
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dissipation, while IitoiUJi — Espm — Wa/{1 — a) becomes much larger than that, the 
more a gets close to unity. Therefore, neglect of W is actually justified when a is close 
to 1, a condition expected to hold in steady-state. Indeed, the expression of eq. H1U|) for 
1 — ckoo gives a correction ^ 1% to for typical cases of interest. 

4. Model application: equations of practical use 

Eq. H1U|) allows to relate directly the model parameter (1 — aoo) to measured quantities, 
in view of the application of the model to real sources. By contrast, a purely numerical 
treatment would require solving the evolution equation for a with given (unknown) ini- 
tial conditions (asynchronism, orbital period, component masses and primary magnetic 
moment), deriving all other quantities as functions of this solution and comparing to 
observations. 

We have used this latter approach to calculate (1 — aoo) in the numerical solution for 
(1 — a): it turns out that, once the system reaches the appropriate regime, the two 
quantities are essentially coincident. Their ratio at sufficiently short orbital periods^ is 
~ 1.001 — 1.002, decreasing with the orbital period down to P '--^ 100 s or even less. 
Therefore, (1 — aoo) represents the actual degree of asynchronism of systems with suf- 
ficiently short orbital period. It can thus be used to obtain physical constraints on the 
two systems directly from measured quantities. 

4.1. The asynchronism evolution timescale 

We have shown that the value a^c (and, accordingly, all steady-state quantities) can be 
determined in full generality from energetic considerations and can be related to measured 
quantities and fundamental system parameters (as described below). The evolution of a 
system in a short-lived, highly asynchronous regime, can be determined only from the 
integration of the evolutionary equation. In this case, indeed, the free-energy reservoir is 
represented by the (unknown) primary's degree of asynchronism. 

In paper I we have studied the evolution equation of a in the linear approximation, with 
constant coefficients. This corresponds to the regime where spin-orbit coupling gives just a 
negligible contribution to the orbital evolution and, at the same time, the synchronization 
timescale Tq is <i; Tq. A system with an initial arbitrary value of a may not be in the linear 
regime: however, the UIM always tends to bring the primary into almost synchronous 
rotation, so that the linear regime will eventually apply. At this point, in order for the 
condition <C Tq to hold, (1 — ctoo) must be sufficiently small or, equivalently, the 
orbital period sufficiently short. As shown in paper I, this regime is formally reached for 
1 — ttoo < 3 X 10^^ or, equivalently, when the orbital period is shorter than the critical 

^ typically ~ 1000 s, see the definition of Pfast in the next section. 
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period: 



0.9 Mq J V 2-5 X 1030 J V 6 X 108 

i?2 \ \^2. 8x10^50 J 



(11) 



.1.7x10V g(i + ,)i(^)Xe)- 

where scaling factors are in c.g.s. units. Pfast measures the period at which, in the nu- 
merical solution, (1 — a) is very well approximated by (1 — ccoo)- 

Strictly speaking, then, if the orbital period is below the Hmiting value Pfast the linear, 
time-independent approximation applies and a reaches ccoo on a time shorter than ~ 0.03 
times the orbital evolutionary timescale. The probability of catching a system during the 
short-lived, asynchronous phase is accordingly very small. Most systems with P < Pfast 
are caught in the long-lived steady-state. 

We finally note that, when a is very close to 1, it can change only by a factor <C a. 
Therefore, the timescale Ta = (x/a can become very long just because of this. However, 
as a passes from 1.1 to 1.001, changing by ~ 10% only, (1 — a) decreases by two orders of 
magnitude and the same holds for E^^''\ while the dissipation rate W changes by four 
orders of magnitude. 

In this situation, the evolution timescale of 1 — a provides more direct information on 
the evolving properties of the system, because this parameter evolves towards zero. We 
thus define, neglecting the minus sign in the last equality: 

1 — a.,~ 



1-a 1-a 
^ ' a(l — a) I at a 



a 



l-a 



(12) 



The difference between this and Tq, is negligible as long as a » 1 and becomes important 
when |1 — a| < 1. In particular, for a = aooj = aoo/(l — Q!oo)to, much longer than 
To- Indeed, a changes only very slightly once it is close to 1. 

On the other hand, T(i_q,^) = To, from which the basic meaning of steady-state is 
recovered clearly: it is the regime in which the system's degree of asynchronism evolves 
on exactly the same timescale of the orbital period. 

For any a ^ aoo the corresponding timescale becomes quickly much shorter than To, 
thus hampering the possibility of observing systems that are too far from steady-state, 
as stated above. 



4.2. Relating measurements to system parameters 

We can now summarize the formulae that allow a direct comparison with the measured 
timing properties and inferred luminosities of the two sources. 

The key parameter (1 — aoo) depends only on the measured quantities d>o and lUo (with 
no assumption on the relative strength of GW and spin-orbit coupling), on the primary 
moment of inertia Ii and on the unknown parameter k, the value of which is primarily 
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determined by /ii. Therefore, based on measurements of Woji^o and W, we can construct 
the quantity a^c as a function of k and use it to estimate ^i. In particular, given the 
definition of a^o and W fea. 1101 and ISI respectively), we have 

W^hujoilJo-, -■ (13) 

The source steady-state luminosity can thus be written as: 

Woo = huJot^oil - aoo) (14) 

The equation for the orbital evolution ^ provides a further relation between the three 
measured quantities, component masses and degree of asynchronism. We re- write it here 
in a form that makes the dependence on measured quantitites and unknown parameters 
explicit: 

W 

Egr - g{u)o){ioo/^o) = 71 r (15) 

(1 - a) 

that becomes, inserting the appropriate expressions for Egr and g{uJo) (cfr. §121), and 
neglecting the small correction to g((^o) due to tidal locking of the secondary: 



o .y3- - (16) 

where X = M^^'^qjil + qY^'^ = J^^/^, A4 being the system's chirp mass. 

With all above equations we can now turn to each of the two systems and assess which 

region of the UIM parameter space applies to each of them. 

5. RX J0806+15 

Israel et al. (2003) measured an on-phase X-ray luminosity (in the range 0.1-2.5 keV) 
Lx = 8 X 10'^^((i/200 pc)^ erg s^^ for this source. These authors suggested that the 
bolometric luminosity might even be dominated by the (unseen) value of the UV flux, 
and reach values 5-6 times higher ( I.e. Lbol 5 X 

^q32 gj.g g 1 c?2^p). The optical flux is 
only 15% pulsed, indicating that most of it might not be associated to the unipolar 
inductor mechanism (possibly the cooling luminosity of the white dwarf plays a role). 
Given these uncertainties and, mainly, the uncertainty in the distance to the source, we 
assume here a luminosity W ~ 10'^^(rf/200 pc)^ erg s^^. In §[3we consider the effect of 
assuming a larger source luminosity for RX J0806-I-15 (and RX J1914-I-24 as well), and 
show that conclusions are affected only weakly. 

Israel et al. (2004) and Strohmayer (2005) obtained independent, phase-coherent tim- 
ing solutions for the orbital period of this source over a ~ 10 yrs baseline, that are 
fully consistent within the errors. The solution reported by Israel et al. (2004) is 
P = -3.67(1) X 10-" and P = 321.53033(2). These give loo 2.23 x lO^^^ rad s'^ 
and uJo — 0.0195 rad s^^. 

In Fig. n (see caption for further details), the dashed hne represents the locus of points 
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in the M2 vs. Mi plane, for which the measured lOq and lOo are consistent with being due 
to GW emission only, i.e. if spin-orbit coupling was absent {a ~ 1). This corresponds to 
a chirp mass M ~ 0.3 Mq. 

Inserting the measured quantities in eq. (|13|l and assuming a reference value* of Ii = 
3 X 10^° g cm^, we obtain: 

1-aoo " 1.3 X 1034 - ^ ''^oo- (17) 

In principle, the source may be in any regime, but its short orbital period strongly 
suggests that it is beyond the limiting period Pfast of eq. Ullfl . Therefore, the most likely 
and natural scenario is that it has already settled into steady-state and a = aoo. 
In this case, we see from eq. (|17f) that (1 — aoo) — 8 x 10^^, Once UIM and spin- 
orbit coupling are introduced, the locus of allowed points in the M2 vs. Mi plane is 
somewhat sensitive to the exact value of a: the solid curve of Fig. ^ was obtained for 
a = = 0.992. Ml must be smaller than 1.1 Mq in order for the secondary not to fill 
its Roche lobe, thus avoiding mass transfer. Note that, if Mi > 1.2 M©, Ii < 10^" g cm^ 



RX J0e06+15 




o.a 0.4 0.6 0.8 1 i.a 

Ml (Me) 

Fig. 1. M2 vs. Ml plot based on the measured timing properties of RX J0806-I-15. The 
dashed curve is the locus expected if orbital decay is driven by GW alone, with no spin- 
orbit coupling. The solid line describes the locus expected if the system is in a steady- 
state, with (1 — a) = (1 — aoo) — 8 x lO"^, TJie horizontal dot-dashed line represents the 
minimum mass for a degenerate secondary not to fill its Roche-lobe at an orbital period 
of 321.5 s. Dotted lines are the loci of fixed mass ratio. 

and assuming W — Woo would imply 1 — aoo > 8 x 10^^ jn this case, the solid curve 
would be shifted downwards, crossing the minimum-mass horizontal line at Mi < 1 Mq. 
In turn this would rule out the high Mi values based on which the curve was obtained. 

Possible values of Ii range in a limited interval: Ji — 3.6 x lO^^g cm^ for AI = 0.2 Mq, 
h = 3.75 x 10^° g cm^ for M = 0.5 Mq, h ^ 3.2 x 10^" g cm^ for M = 0.8 Mq, h = 2.4 x 10^'' 
g cm^ if Ml = 1 Mq and h ^ W''° g cm^ for Mi = 1.3 Mq. 
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Therefore, we conclude that, if RX J0806+15 is interpreted as being in the UIM steady- 
state, M < 1 Mq is necessarily required. 

From (1 — Uoo) = 8 x lO^'^ we obtain, by use of eq. k ^ 7.7 x 10^^ (c.g.s.): from 
this, component masses and primary magnetic moment can be constrained (see Fig.^). 
Indeed, since k = k{fii, Mi, q;a) (see its definition in eq.|3|), we have a further constraint 
on it, because Mi and q must lie along the solid curve of Fig.^ Given the value of ct, /ii 
is obtained for each point along the solid curve. 

According to Wu et al. (2002), and references therein, the electrical conductivity of a 
white dwarf atmosphere with a temperature ~ 10^ K has a value ^ lO^'^ ~ 10^^ (e.s.u.). 
We take ct = 3 x 10^^ (e.s.u.), so that the whole range above is within a factor 3. 
The values of fj,i obtained in this way, and the corresponding field at the primary's 
surface, are plotted in Fig. |2 from which /ii ~ a few xlO'^*' G cm"^ results, somewhat 
sensitive to the primary mass. 

We note further that the solid curve of Fig. ^ is not strictly a locus of constant chirp 
mass. However, it is nearly so: X ~ (3.4 4.5) x lO^"' g^^^ along it, which implies 
M ~ (0.26 ^ 0.31) Mq. More importantly, E^, ~ (1.1 ~ 1.9) x 10^^ erg s^^ and, since 
W/{1 - aoo) = e'i''''^ ~ 1.25 X 10^4 erg s'^, we have E^, ~ {9 ~ 15) Therefore, 
orbital spin-up is essentially driven by GW alone; indeed, the dashed and solid curves are 
very close on the M2 vs. Mi plane. Summarizing, we have shown that RX J0806+15 has 




0.4 0.5 0.6 0.7 O.B 0.9 1 1.1 

Fig. 2. The value of the primary magnetic moment /ii, and the corresponding surface 
B-field, as a function of the primary mass Mi, for (1 — a) = (1 — aoo) = 8 x 10^'^. 

observational properties that can be well interpreted in the UIM framework, assuming it 
is in steady-state. This requires the primary to be moderately magnetic, with /ii ^ 10^*^ 
G cm'^ and with an almost synchronous spin, slightly slower than the orbital motion (the 
difference being less than 1%). The small degree of asynchronism is consistent with 
the fact that no modulation of pulse arrival times at the beat period Ub = tOo — ^1 is 
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observed, since its small amplitude would likely be below the intrinsic timing noise of 
the source (cfr. Barros et al. 2004). The expected value of the magnetic moment can also 
be tested by future observations, for example through studies of polarized emission, at 
optical and/or radio wavelenghts (Willes & Wu 2004). 

6. RX J1914+24 

The luminosity of this source has been subject to many revisions over the last few years. 
Wu et al. (2002) refer to earlier ASCA measurements that, for a distance of 200-500 pc, 
corresponded to a luminosity in the range (4 X 10^3 ^ 2.5 X 10^4) erg s^^. Ramsay et al. 
(2005), based on more recent XMM-Newton observations and a standard blackbody fit 
to the X-ray spectrum, derived an X-ray luminosity of ~ lO^^^^fepc where dkpc is 

the distance in kpc. The larger distance of ~ 1 kpc was based on a work by Steeghs et 
al. (2006). Still more recent XMM-Newton observtions have been analyzed by Ramsay 
et al. (2006): these authors find that an optically thin, thermal emission spectrum, with 
an edge at 0.83 keV (attributed to O VIII) gives a significantly better fit to the data 
than the previously used blackbody model. Compared to a blackbody, the optically thin 
thermal plasma model clearly leads to a source bolometric luminosity much lower than 
previously derived, Lboi — 10^'^ d^^^ erg s^^. Ramsay et al. (2006) also note that the 
determination of a 1 kpc distance is not free of uncertainties and that a much smaller 
distance 200 pc) is still possible: based on the latter, the luminosity would become as 
small as 3 X 10"^^ erg s^"'^. 

Given these large discrepancies, interpretation of this source's properties is clearly am- 
biguous and quite dependent on which assumptions are made. We refer here to the more 
recent assessment by Ramsay et al. (2006) of a luminosity L = 10'^^ erg s^^ for a 1 
kpc distance, but want to stress mainly the general characteristics of our procedure. The 
numerical estimates that we obtain are, at this stage, necessarily illustrative and, in §[3 
we show how derived quantities change when different luminosities are assumed. 
Ramsay et al. (2006) also find possible evidence, at least in a few observations, of two 
secondary peaks in power spectra. These are very close to (Av ~ 5 x 10~^ Hz) and sym- 
metrically distributed around the strongest peak at ^ 1.76 x lO^'^ Hz, the latter usually 
interpreted as the orbital period of the source. We will briefly comment on these possible 
features at the end of § 16.31 

Ramsay et al. (2005), confirmed by Ramsay et al. (2006), also give the most up to date 
measurement of the source's rate of orbital evolution, P — 3.2(10) x 10^^^, that converts 
to ujo — 6.2 X 10^^'' rad s^^ for the orbital period of 569 s, or cJq — 0.011 rad s^^. We 
warn that these parameters were not obtained, like in the case of RX J0806+15, through 
a phase-coherent analysis, therefore they may be subject to larger uncertainties and fu- 
ture revision. 
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Application of the scheme used for RX J0806+15 to this source is not as straightforward, 
since its properties indicate that it is not in a steady-state and, hkely, not even in the 
regime where spin-orbit couphng is neghgible in the orbital evolution. 
First of all, its inferred luminosity seems inconsistent with being associated to steady- 
state as, with the measured values of Uq and Uo, eq. H13|l gives^ 

(1 - a? ^ 10'' ,2 
1 - aoo 2 X 1032 '^'^■f^ 

- 5dlpc . (18) 

which has important implications. First of all, given ujo and uJo, the system steady-state 
luminosity should be ^ 2 x 10^^ erg s^^ fea. 114(1 . hardly consistent with measurements. 
Indeed, even if a large value of (1 — ctoo) — 10~^ was assumed, ~ 2 x 10^^ erg s~^, 
less than the smallest possible luminosity reported by Ramsay et al. (2006). A relatively 
high ratio between the actual asynchronism parameter and its steady-state value appears 
unavoidable. Indeed, from the above equation: 

|l-a| -2.2(1 -aoo)'/' (19) 

For example, if this system had the same k and Ii as RX J0806-I-15, it would have a 2.5 
times smaller value of (1 — aoo) — 3.3 x 10"^^ from which |1 — a| ~ 0.13. 



6.1. The case for a > 1 

The low rate of orbital shrinking implied by the measured lUo — 6.2 x 10^^^ s^^ (Ramsay 
et al. 2005) suggests a very low mass system, with Mi < 0.3 Mq; in this case the 
secondary would not fill its Roche- lobe only for q > 0.26 (sec Fig.lJand Nelemans 2004 
for a discussion of this point). In the UIM framework, this result gives important hints 
as to the primary spin. 

If, as shown above, the system is far from steady-state and has a < 1, the locus of 
allowed masses in the M2 vs. Mi plane is shifted downwards with respect to the dashed 
line of Fig.|3|by a significant amount. This would exacerbate the problem of the very low 
component masses. In particular because, a being quite different from 1, the downward 
shift of the curve in the M2 vs. Mi plane would accordingly be large. 
This argument can be made more cogent by using the measured value of > 0; a strong 
constraint on |1 — a| results, because the high inferred luminosity W ^ lO^^ erg s~^ 
would imply a high value of E^^^^'K If a < 1, this latter term sums to the GW torque: the 
resulting orbital evolution would thus be faster than if it were driven by GW alone. In 
fact, the smallest possible value of E^^^^'^ for a < 1 is obtained with a = 0, from which 

Therefore, if a < 1, an absolute minimum to the system's rate of orbital shrinking would 
^ again assuming /i = 3 x 10^" g cm^ 
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result from the above minimum rate of orbital energy extraction operated by spin-orbit 
coupling. This minimum ibo turns out to be very close to the measured one, so that it 
would imply unplausibly small component masses in order for E^^ to be sufficiently small. 
This can be shown by comparing the corresponding chirp masses. The dashed curve in 
Fig. m corresponds to a value of X ~ 1.07 x 10^"' g^^'^, or a chirp mass Ai ~ 0.13 Mq. 
On the other hand, an equivalent GW luminosity of lO'^'^ erg s^^ would correspond to 
X ~ 8.5 X 10^^ g^/^ or an equivalent chirp mass M ~ 0.12 M©. 

Therefore, the slow measured rate of orbital evolution of this source essentially rules out 
the case a < 1 in the UIM framework discussed here: indeed, no point in the M2 vs. 
Ml plane can be found for a < 1 if a non-Roche-lobe-filling secondary is required (as 
it is in the UIM). We note that the absolute minimum is obtained assuming a = 0, a 
very unlikely condition. Summarizing, Fig. |21 and the above argument strongly lead to 
consider the case where the primary spin is faster than the orbital motion. This would 
indeed offer a great advantage in interpreting this source: in the case a > 1, spin-orbit 
coupling has an opposite sign with respect to the GW torque. The two mechanisms have 
opposite effects on the orbital evolution and a slow measured ujo can result from the two 
opposite torques partially cancelling each other. 

We want to stress further the above arguments, as there appears to be some confusion 
in the literature. Marsh & Nelemans (2005) have recently claimed that, in order for the 
UIM to account for the properties of RX J1914+24, a > 10^ would be required, implying 
the primary is spinning nearly at break-up. 

These authors arrive at their conclusion first deriving an estimate of the GW luminosity 
from the measured ujo and lOq and then requiring that W/{l—a) be less than the estimated 
-Egr. Given the high source luminosity (they assume 10^^ erg s~^) and small uJo, a large 
value of |1 — a| is required in order for that condition to hold. Note that, assuming W~ 
2^033 gj.g g-i^ their argument gives |1 — a| > 1, not so extreme indeed. Most importantly, 
however, this argument relies on the assumption that the measured timing parameters of 
this source give a reliable estimate of its actual GW luminosity. This assumption holds 
only if spin-orbit coupling is negligible; we have shown above that this needs not be the 
case and, most likely, is not the case in this source. The argument by Marsh & Nelemans 
might thus be reversed and used as evidence that a non-negligible spin-orbit coupling is 
necessary in applications of the UIM to this source. 

Ramsay et al. (2006) suggest that the much smaller luminosity of lO'^'^ erg s^^ they 
derive would weaken the objection by Marsh & Nelemans, since it would require an 
asynchronism (1 — a) ^ 10~^ only (or even less). It is easy to see that this is not true, 
because with W= 10^^ ergs^^ and (1 — a) < 10^^ one would obtain = W/{l-a) > 

10'^^ erg s^^, which is greater than any plausible GW emission from this system, given its 
measured uJo- Rather, |1 — a| cannot be too small without making the model inconsistent 
with the observed slow orbital spin-up. 
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The conclusion of this section can only be that, in a self-consistently constructed UIM 
framework, the GW luminosity of this source has to be larger than implied by uJo, which 
implies a > 1, otherwise the model hardly applies. 



6.2. Constraining the asynchronous system 

We try now to give a quantitative description of the possible scenario introduced above. 
The goal is to find constraints on system parameters in order to match the measured lu- 
minosity and uJo/lUo and meet the requirement that the resulting regime has a sufficiently 
long lifetime. Given that it is not possible to determine all system parameters uniquely, 
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Fig. 3. M2 vs. Ml plot based on measured timing properties of RX J1914+24. The dot- 
dashed line corresponds to the minimum mass for a degenerate secondary not to fill its 
Roche-lobe. The dashed curve represents the locus expected if orbital decay was driven 
by GW alone, with no spin-orbit, coupling. It is seen that this curve is consistent with a 
detached system only if both masses are extremely low. The solid lines describe the loci 
expected if spin-orbit coupling for the primary star (the secondary spin is tidally locked 
to the orbit) is present and gives a negative contribution to lUo- The four curves are 
obtained assuming W = 10'^^ erg s~^ and four different values of a = 1.025, 1.05, 1.075, 
1.1, respectively, from top to bottom, as reported in the plot. 



the scheme that we have used is as follows: given a value of a, the measured Woi^o and 
W allow us to determine a from eq. I@J. For each value of Mi eq. re- written in the 
form reported in the Appendix, can thus be used to determine the value of M2 (or q) 
that is compatible with lUo and Wq. This yields the solid curves of Fig. |21 
As these curves show, the larger is a and the smaller the upward shift of the corre- 
sponding locus, which may seem surprising. However, these curves are obtained at fixed 
luminosity W and luo- Recalling that (l/a) gives the efficiency of angular momentum 
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transfer in systems with a > 1 (cfr. a higher a at a given luminosity imphes that 

less energy and angular momentum are being transferred to the orbit. Accordingly, GWs 
must be emitted at a smaller rate for a given ujo. 

The values of a in Fig. |31 were chosen arbitrarily and are just illustrative: note that the 
resulting curves are quite similar to those obtained for RX J0806+15. Given a, one can 
also estimate k from the definiton of W (eq.OJ. On the other hand, the information given 
by the curves of Fig. yields determination of all quantities contained in k, apart from 
/xi. Therefore, assuming ct = 3 x 10^^ (e.s.u.) as in the previous case, we can determine 
the value of /ii along each of the four curves of Fig. Derived values are plotted in Fig. 
^ We see that /ii varies by just a factor 2-3 as a function of Mi along each curve, while 
it depends somewhat more strongly on the assumed value of a. 

The resulting numbers correspond in any case to a moderately magnetic primary, indeed 
quite similar to what obtained for RX J0806+15 (the surface field B is in the 10^ lO'^ 
G range, depending on the primary mass). 



RX J1914+24 
W = 10^^ erg 




5 0.4 0.6 0.8 1 1.2 
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Fig. 4. Ill as a function of the primary mass for the same values of a used previously 
and reported on the corresponding curves. Given a value of a and the estimated lumi- 
nosity W ^ 10^'^ erg s~^, a is calculated from eq. Q as a function of Mi. q as well is 
obtained as a function of Mi from Fig. 13 All system parameters contained in k are thus 
determined, apart from /ii. The estimated luminosity W and the assumed value of a 
give the corresponding value of k: combining all the information, fii is self-consistently 
determined. 



6.3. Lifetime of the asynchronous state 

Finally, having argued that this system can be explained in the UIM framework by a 
relatively highly asynchronous and moderately magnetic system, we have to deal with 
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the timescale problem. 

While RX J0806+15 is fully consistent with being in steady-state, so that it has no 
timescale problem, we have demonstrated that RX J1914+24 cannot be in the same 
regime. Moreover, given the combination of a high luminosity and slow orbital evolution, 
we have found arguments in favour of its primary spinning faster than the orbital motion. 
It must be fast enough, indeed, that the orbital evolution is significantly slowed down 
by the angular momentum transferred to the orbit by spin-orbit coupling. Therefore, the 
system is in the fully non-linear regime where spin-orbit coupling is non negligible and 
the evolution equation for lOq must be solved explicitly. 

We can estimate the synchronization timescale starting from the evolution equation 
(0J. With the measured values of W , lOq and Wq, can be calculated as a function of Ii 
(and, thus, of Mi) given a particular value of a. In Fig. [S] we have reported the results 
obtained for the same four values of a assumed in the previous section. 
It is apparent that the larger a, the longer the timescale, for a given value of Mi. This 
is because the system luminosity is a fixed (measured) quantity: therefore the larger 
the asynchronism, the smaller the efficiency of spin-orbit coupling and, accordingly, the 
larger Tq. The resulting timescales range from a few xlO^ yrs to a few xlO^ yrs, while 
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Fig. 5. The evolution timescale Tq as a function of the primary mass for the same values of 
a used previously, reported on the corresponding curves. Given the estimated luminosity 
W ~ 10'^^ erg s~^ and a value of a, Ta is calculated from eq. as a function of h and, 
hence. Mi. 

To 6 X 10^ yrs for this source. Although may still seem too short, we stress that 
in this scenario the long orbital evolutionary timescale results largely from the effect of 
spin-orbit coupling. The latter decreases strongly on the timescale Tq, after which Tq is 
essentially determined by GW emission only. We denote the latter as r^^. Therefore, 
the system will actually spend the total time at the nearly constant orbital period. 
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its evolution being slowed down by just a fraction t^/t^^ of this time. In Fig.[Slwe have 
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Fig. 6. The ratio r^^^/ra calculated for the four values of a = 1.025,1.05,1.075,1.1. 

is obtained from the value of X along each of the four curves of Fig. |31 according 
to eq. (|20|l . while the corresponding value of Tq, are those plotted in Fig.|S| 



plotted the ratio (t^'^/tq) for the four cases studied above (Figs.|3IElandEJl- The values 
were calculated along the solid curves of Fig.|21 that define the value of X for each 
(Mi,M2). Given X, both the source GW luminosity E^-, and t^"^ can be calculated. In 
particular, since gi^Oo) is a function of X itself, we have 

GW gg ^5/3 



-) = (rf = . (20) 

LOo } b & 

According to the above reasoning, the expected GW luminosity of this source is in the 
range (4.6 1.4) x 10'^'* erg s~^. The corresponding ratios E^^/E^l^''^ 1-1^^ 1.21, 1.29 
and 1.4, respectively, for a ~ 1.025, 1.05, 1.075 and 1.1. 

Therefore, RX J1914+24 and RX J0806+15 may even have similar component masses 
and magnetic moments. The significant observational differences would essentially be due 
to the two systems being caught in different evolutionary stages. RX J1914+24 would be 
in a luminous, transient phase that preceeds its settling into the dimmer steady-state, a 
regime already reached by the shorter period RX J0806+15. Although the more luminous 
phase is transient, its lifetime can be as long as ~ 10^ yrs, one or two orders of magnitude 
longer than previously estimated. 

We comment here on the two symmetric peaks, around the main one, found in a few 
power spectra by Ramsay et al. (2006). 

In the UIM framework, given a primary spin ui that differs from ujo, some signature of 
both periodicities may be expected, indeed. The hot spot position on the primary surface 
is modulated at the orbital period (the spot always follows the position of the secondary 
star) and at the beat period ujb = oji — oJq- The latter because the spot position is deter- 
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mined by the position of the primary's magnetic pole as well, which itself rotates at lui. If 
the field simmetry axis were perfectly aligned to the spin (and orbital) axis, this second 
modulation would disappear but, in general, some inclination is to be expected. Barros 
et al. (2004) have shown that, in the context of the UIM, the modulation at ujb in RX 
J1914+24 is small, possibly negligible, if the inclination angle f3 < (15 — 19)°. However, 
it would be a natural expectation of the model to have two main periodicities in power 
spectra, namely Wo and ujb: when uji ujo (thus a ^ 1) the frequency Ub Wo- 
This is highly speculative at this stage and proper quantitative studies of this phe- 
nomenon, both on the observational and theoretical side, are necessary before these 
features can be used to put quantitative constraints on models. These are beyond the 
scope of this paper and are currently in progress. What we want to stress here is that 
a second periodicity in power spectra is not a finding against the UIM. The features 
reported by Ramsay et al. (2006) are not in constrast with qualitative expectations in 
the UIM. 

7. Dependence of the results on the assumed luminosity 

In this section we briefly discuss how our results change when different source luminosities 
are assumed, reflecting the large uncertainties in their distances and bolometric emission. 
We have repeated the same reasoning presented in the previous sections, simply changing 
the adopted value of W to a significantly larger value. 

7.1. RXJ0806+15 

We have discussed in fj^the uncertainty in the bolometric luminosity of this source, which 
could possibly be up to a factor 5-6 larger than the X-ray emission reported by Israel et 
al (2003). We assume here W =5 x 10'^^ erg s~^ and see how our conclusions change. 
For this luminosity, a six times larger degree of asynchronism than previously assumed 
would result, 1 — ctoo — 3.8 x 10^^, still a quite small number. Given this value, the 
corresponding locus in the M2 vs. Mi plane would correspond to values of X and, ac- 
cordingly, GW luminosities smaller by just a few percent. The corresponding primary 
magnetic moment and implied surface field are themselves decreased by the same factor 
of six by which W increases. All these changes, although non-negligible on quantitative 
grounds, do not affect significantly the overall picture presented above. 

7.2. RX J1914+24 

The case of RX J1914+24 is more difficult: the luminosity assumed in § lies at the 
centre of a very wide range of proposed values, from the smallest W ~ 3 x 10^^ erg s~^ 
(Ramsay et al. 2006) to the largest W ~ 10'^^ erg s^^ (Ramsay et al. 2005). We show 
and briefly discuss here the same figures presented in the previous section, obtained for 
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Fig. 7. The loci in the M2 vs. Mi plane obtained assuming a luminosity of W = 3 x 10^"'^ 
erg/s and a = 1.005 (solid line) or a very high luminosity of ^ 10'^^ erg/s and a — A 
(dashed line). As noted in the previous section, curves shift upwards for smaller values 
of a and downwards for increasing values of the asynchronism. 

just the two extreme luminosities above. Instead of exploring a range of values of a, we 
are exploring the range of values of W , over approximately four orders of magnitude. 
Therefore just one curve, is reported for each of the two cases. In all figures below the 
solid lines correspond to the low luminosity W=: 3 x lO'^^ erg s^^, while the dashed 
curves represent the opposite case with W= lO'^^ erg s~^. From Fig. Qwe see that the 
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Fig. 8. The ratio t^^ /tc, calculated for W= 3 x lO-^^ erg/s and a = 1.005 (solid line) 
and for W= 10"^^ erg/s and a — A (dashed line), is obtained from the value of X 

along the corresponding curve of Fig. according to eq. (|20|l . while the corresponding 
value of Ta are calculated as described in the caption to Fig. O 



high luminosity case would require a significant degree of asynchronism, but still of the 
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order of a few, not the ^ 100 value previously suggested in the literature. On the other 
hand, the very dim case can be well described by a low asynchronism, even lower than 
the steady-state asynchronism found for RX J0806+15. However, we have shown that 
1 ~ oioa — 0.15 would be implied for this source if 3 x 10"^^ erg s^^ — Woo- Rather, 
a = 1.005 gives {1 - Uoo) ^ 1.7 x IQ-* (eq. dl. 

The dimmer system would have a long lifetime, Tq 3 x 10^ yrs, which is ^ one tenth 
to one twentieth of the orbital evolutionary timescale ■ From the solid curve of Fig. 
[7|we can estimate a GW luminosity E^^ ~ 10'^'* erg s^^ and a ratio E^^/ E^f^'^^^ ^ 1.6. On 
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Fig. 9. The expected value of the primar's magnetic moment /ii as a function of the 
primary mass (Mi) for the same two extreme luminosities of the above figures, indicated 
in the figure. Despite the non negligible difference between the two curves, these still 
span the range of medium-to-weak magnetization already found previously. 

the other hand, the very luminous system has a shorter Tq, ~ a few x 10** yrs, but can still 
have a lifetime ^ one tenth of the orbital evolutionary timescale ■ The latter is itself 
significantly shorter than ojo/wo, because large component masses would be required and, 
correspondingly, very strong GW emission expected. For the case with a = 4 we estimate 
Eg^ ~ 3.6 X lO"^"* erg s^^, with a ratio ^ 1.25 to the rate of work W/{1 — a). 
Finally, Fig. |3 shows the expected primary magnetic moments in both cases. 

8. Conclusions 

In the present work we have clarified some key aspects of spin-orbit coupling in the 
Unipolar Inductor Model coupled to GW and applied the model to the two candidate 
ultrashort period binaries RX J0806-I-15 and RX J1914+24. We derive constraints on 
system parameters in order for the model to account for the inferred luminosities and 
measured ujo and ujo of both sources. The UIM accounts in a natural way for their 
luminosities, measured orbital periods and rates of orbital spin-up, provided they have 
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moderately magnetic primaries, say fii ~ several x 10^^ to few x 10^° G cm'^. Concerning 
each of the sources, we conclude that: 

— In the case of RX J0806+15 the observations are explained in the framework of 
the steady-state solution. The degree of asynchronism is determined by a balance 
between the rate of current dissipation and the rate at which energy is fed to the 
electrical circuit by orbital evolution, while the source luminosity does not decrease 
at all. The orbital evolution is driven mostly by GW-emission and the system's timing 
properties should differ only slightly from those of two point masses evolving under 
GW emission. The modest primary magnetic moment represents a key requirement 
of this interpretation and its measurement would be a crucial test as to its viability. 

— RX J1914+24, on the other hand, cannot have reached the steady-state solution yet, 
as its present rate of electric energy dissipation is significantly higher than what GW 
emission alone could feed to the circuit. 

In addition to requiring a moderately magnetic primary (again ^ 10^° G cm^), we 
conclude that the primary spin frequency must be somewhat larger than the orbital 
frequency. For W~ 10^^ erg s~^ the parameter (1 — a) can be in the range (0.1 -j- 0.02) 
(Pspin between 517 and 558 s). Although the source is not consistent with being in a 
steady-state we have shown that, for a in the above range, the asynchronous, luminous 
and transient phase can last between a few xlO* yrs up to 3 x 10^ yrs. These are 
between 2% and 10% of the orbital evolution timescale. Therefore, there is a sizeable 
probability that one in two of the ultrashort period binaries that have so far been 
found is in this evolutionary stage. 

The measured value of Wo does not reflect the actual rate of GW emission from this 
system. It rather results from the sum of GWs and spin-orbit coupling, the latter 
partially cancelling the effect of the former. For the range of values of a that we 
consider more likely, a ~ (1.025 1.1), the actual GW luminosity should be in the 

^ lO'^'* erg s~^ range. 

Values of a between 1.005 and ^ 4 5 can explain all possible source luminosities in 
the range (3 x lO''^ lO'*"'') erg s^^. In particular, a ~ of a few can self-consistcntly 
explain a luminosity as high as 10^^ erg s~^, extremely less then previously claimed 
in the literature. 

In this work we were mainly concerned with finding system parameters for which the 
model could be made consistent with the energetic requirements and timing properties 
of the candidate sources. Further work is needed to address issues such as the duration 
and shape of the X-ray pulsations. Predictions that can be checked through observations 
to be carried out in the next few years are essentially of three kinds: 

— independent measurements of the primary magnetic moment y^i, which should result 
to be smaller than several x 10^° G cm^ in order for the model to account for the two 
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sources. Although several details of the model should still be taken into account, it 
seems difficult that the overall energetic requirements could be changed by orders of 
magnitude. The expectation of low magnetic moments seems therefore quite robust. 

— extending the temporal baseline over which timing studies are carried out may help de- 
tecting the second derivative of ojo in both systems, but particularly in RX J19f 4+24. 
Such measurements could provide an indipendcnt estimate of system masses and/or 
confirm a significant contribution of UI to the orbital evolution. 

— a detailed study of the duration and shape of the X-ray pulsations, as well as of power 
spectra, would be required, in order to put constraints on geometric properties of the 
systems, possibly related to the shape and orientation of magnetic field lines near the 
primary component. 

Account for hot spot shapes different from a simple filled spot should be taken since 
two detached, arc-like hot regions are expected instead (see Wu et al. 2002). 

Concerning the longer term future, we emphasize that LISA is expected to measure the 
GW emission of very short period binaries such as those discussed here, therefore yielding 
an independent check of the model. In particular, RX J1914-f24 is expected to be much 
more luminous in GWs than one would infer based on its measured lOq and lOq. 



9. Appendix 

We derive here the expression used to obtain the solid curves of Fig. |2| In the case of 
Fig. n a was assumed to be constant, as this describes the relevant condition to RX 
J0806-I-15. On the other hand, in studying RX J1914-I-24 we were led to deal with the 
more general case in which a has a non negligibile time derivative. 
Let g{ijJo) be defined as: 
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where the second term in square brackets is denoted by f{u!o)- Then, from eq. ||2Jl, 
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and from this, re-writing the pure-GW period derivative as {ujo/ujo)g — Egrf{uJo)/gi^o)- 
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that leads to the following expression: 



{uJo/uJo). 



gi^o) 



f{uJo) {uJo/LOo)g 



(24) 



24 S.DaU'Osso et al.: Astrophysical UI powered by GW emission 

Now substitute the expressions for g{u!o) and f{i^o), remembering that Ii = (2/5)Mii?J 
and that, from Kepler's third law, ujo^^^ = [GMi{l + q)]^^^'^a'^. After some algebra it 
obtains: 
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The solid curves of Fig. |31 were obtained just using eq. (|25|l . This takes account of the 
primary asynchronism and its rate of variation as well, while assumes the secondary to 
be perfectly synchronous (tidally locked) at all times. 
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